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Abstract 

Background: Windthrow, the uprooting of trees during stornns associated with strong winds, is a well-established 
cause of nnortality in temperate regions of the world, often with large ecological consequences. However, this 
phenomenon has received little attention within arid regions and is not well documented in southern Africa. Slow 
rates of post-disturbance recovery and projected increases in extreme weather events in arid areas mean that 
windthrow could be more common and have bigger impacts on these ecosystems in the future. This is of concern 
due to slow rates of post-disturbance recovery in arid systems and projected increases in extreme weather events 
in these areas. This study investigated the spatial pattern, magnitude and likely causes of windthrown mortality in 
relation to other forms of mortality in Aloe dichotoma, an iconic arid-adapted arborescent succulent and southern 
Africa climate change indicator species. 

Results: We found that windthrown mortality was greatest within the equatorward summer rainfall zone (SRZ) of 
its distribution (mean = 31%, n = 1 1), and was derived almost exclusively from the larger adult age class. A logistic 
modelling exercise indicated that windthrown mortality was strongly associated with greater amounts of warm 
season (summer) rainfall in the SRZ, higher wind speeds, and leptosols. A statistically significant interaction term 
between higher summer rainfall and wind speeds further increased the odds of being windthrown. While these 
results would benefit from improvements in the resolution of wind and substrate data, they do support the 
hypothesised mechanism for windthrow in A dichotoma. This involves powerful storm gusts associated with either 
the current or subsequent rainfall event, heavy convective rainfall, and an associated increase in soil malleability. 
Shallow rooting depths in gravel-rich soils and an inflexible, top-heavy canopy structure make individuals especially 
prone to windthrown mortality during storms. 

Conclusions: Results highlight the importance of this previously unrecognised form of mortality in A dichotoma, 
especially since it seems to disproportionately affect reproductively mature adult individuals in an infrequently 
recruiting species. Smaller, more geographically isolated and adult dominated populations in the summer rainfall 
zone are likely to be more vulnerable to localised extinction due to windthrow events. 

Keywords: Aloe dichotoma, Arid environment, Convective rainfall, Indicator species, Mortality, Southern Africa, 
Wind speed, Windthrow 



Background 

Windthrow, which refers to trees being broken or 
uprooted by the wind, is a well established cause of cata- 
strophic mortality in plantations and natural forests, 
often with large economic and ecological consequences 
[1,2]. Analyses of the impacts of windthrown mortality 
have typically been restricted to tropical forests and 



* Correspondence: timm.hoffman(5)uct.ac.za 

Vlant Conservation Unit, Department of Biological Sciences, University of 
Cape Town, Cape Town, South Africa 

Full list of author information is available at the end of the article 

(3 BioMed Central 



northern hemisphere temperate forests e.g. [3-7]. Very 
few studies have documented the ecological effects of 
windthrown mortality on species in more open, arid en- 
vironments (but see [8] for a savanna example with a so- 
cial focus). 

Understanding the effects of windthrow in arid ecosys- 
tems is important for two reasons. Firstly, plant commu- 
nity recovery from disturbance, such as windthrown 
mortality, is typically more prolonged in arid systems 
due to resource limitations such as a lack of moisture 
and, therefore, access to nutrients. Resource limitations 
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reduce growth rates and extend the intervals during 
which phases of reproduction and estabUshment are 
likely to occur [9,10]. Added anthropogenic pressure in 
arid regions come in the form of landscape fragmenta- 
tion, heavy grazing pressure, and water abstraction, 
which combine to further retard rates of recovery [11]. 

Secondly, global climate change models predict an in- 
crease in the intensity and a decrease in the return inter- 
val of extreme weather events [12]. Notwithstanding 
evidence for a counter- intuitive decline in global average 
windrun [13], an increased intensity of storm events is 
likely to be accompanied by increased precipitation and 
localised wind speeds in certain areas e.g. [14], with im- 
plications for windthrown mortality, as has been postu- 
lated for European forest stands [15,16]. 

The current study investigated the role of windthrown 
mortality in Aloe dichotoma, a widespread and iconic 
succulent tree species distributed in the arid and semi- 
arid winter and summer rainfall regions of Namibia and 
South Africa. Previous research has emphasized the role 
that increasing temperature, drought, and biotic factors 
such as herbivory and disease have played on the high 
levels of mortality observed in many populations of 
this species [17-19]. Despite observational evidence that 
windthrow contributes a large fraction to total mortality 
in some A. dichotoma populations, this phenomenon 
has not been investigated. Indeed, an explicit analysis of 
windthrown mortality does not exist for any southern 
African species, despite the fact that conditions likely to 
promote windthrow events (i.e. convective storm events 
with associated high wind speeds) are ubiquitous in the 
summer rainfall zone (SRZ) which dominates southern 
Africa e.g. [20]. In addition, southern Africa appears to 
be an area where extreme weather events are predicted 
to increase in frequency [21,22]. 

Studies in both temperate and tropical settings suggest 
that taller, larger individuals with greater canopy surface 
areas are likely to be more susceptible to windthrown 
mortality [23,24]. This is due chiefly to greater drag and 
reduced stem flexibility associated with large individuals 
[25,26], although other architectural factors such as 
reduced canopy density, increased branch flexibility, and 
a more aerodynamic leaf-shape can have a mitigatory 
effect [27]. This may be cause for concern when consid- 
ering A, dichotoma^ as the species has a characteristically 
rigid stem, canopy branch and leaf architecture. In 
addition, smaller and more climatically marginal popula- 
tions, which are often dominated by large, mature adult 
individuals, may be especially at risk. Adult individuals 
typically bear a disproportionate amount of the repro- 
ductive potential in a population and windthrow events 
which increase adult mortality can therefore have a det- 
rimental effect on long term population persistence, es- 
pecially in areas where recruitment is infrequent. 



Soil water content, as well as rooting structure, depth, 
and substrate, are important additional determinants of 
windthrow susceptibility. An increase in soil water con- 
tent, for example, reduces soil resistance to slippage 
[28], while an increase in the adhesive ability of soil 
particles to water presumably prolongs the period of 
vulnerability to windthrow [29]. Generally, thicker root 
plates (i.e. the vertical rooting extent) which tap deeper 
water sources in well drained soils provide a greater re- 
sistive force against shearing [30]. The less horizontally 
extensive and shallower root plates typical of adventi- 
tious root systems associated with monocotyledonous 
plants such as A, dichotoma may, therefore, be more 
susceptible to shearing at the root-soil interface, par- 
ticularly if the soil at the base of the root plate (and 
below this level) is saturated. This phenomenon has 
been well demonstrated for the root anchorage of Cha- 
maecyparis obtusa (Hinoki), a cypress species native to 
central Japan [28]. 

Given the general lack of knowledge pertaining to this 
phenomenon in semi-arid southern African, the current 
study provides baseline information on the spatial pat- 
tern and demographic characteristics of windthrown 
mortality in A, dichotoma. It also establishes the per- 
centage of mortality resulting from windthrow, rather 
than from other possible causes, such as physiological 
and climatological drought e.g. [17,19,31]. Lastly, we re- 
late the pattern of windthrown mortality to potential 
causative factors such as tree size, rainfall characteristics, 
wind speed, geology and lithology. In doing so, we iden- 
tify areas where populations might be more vulnerable 
to windthrow, and therefore require more active man- 
agement to ensure their continued survival. 

Methods 

Data collection and classification 

Fourteen A, dichotoma populations were sampled over 
two consecutive winter seasons in Namibia and the 
Northern Cape in South Africa during 2008 and 2009 
(Table 1). As far as possible, populations were chosen at 
roughly equal intervals spanning the latitudinal extent 
and thus the climatic range of the species' distribution. 
Populations with individuals on all four aspects were 
also selected so as to capture the full range in local con- 
ditions due to micro-topographical effects. Where pos- 
sible, approximately 240 (60 per aspect) individuals 
(both living and dead) were recorded along altitudinal 
transects within a population. Conspicuous populations 
along busy transport routes easily accessed from the 
roadside were avoided due to the likelihood of prior re- 
moval of juveniles or (potentially windthrown) dead 
skeletons for firewood or building purposes. 

GPS co-ordinates and altitude were recorded for every 
individual and these were used to extract average climatic 
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Table 1 Numbers of sampled Aloe dichotoma Individuals per population and basic climate characteristics 



Population name 


Total number of 
individuals counted 


Ave. altitude (m) 
(±stdev) 


Ave. rainfall (mm) 
(±stdev) 


Ave. temperature 
rC) (+stdev) 


Brandberg 


177 


1666 (±277) 


204.3 (±41.8) 


17.3 (±1.4) 


Omaruru River 


185 


734 (±25) 


100.8 (±4.2) 


21.8 (±0.1) 


Spitzkoppe 


162 


1231 (±86) 


135.1 (±9.9) 


20.1 (±0.3) 


Tinkas River 


258 


699 (±25) 


82.2 (±2.4) 


20.6 (±0) 


Remhoogte 


211 


1411 (±89) 


1 894 (±9.4) 


16.9 (±0.5) 


Hauchabfontein 


290 


1226 (±55) 


1 24.3 (±0.5) 


1 6.6 (±0) 


Gorab 


Z/H 




1 OQ A (-\-Q 


\ J.J [±vJ.D) 


Namtib 


293 


1265 (±59) 


102.8 (±2.9) 


154 (±0.1) 


Carolinaliof 


222 


1273 (±12) 


134 (±1.5) 


17.2 (±0.1) 


Kliplioek 


293 


1262 (±36) 


101.5 (±4.4) 


15.5 (±0.2) 


Grunau 


241 


913 (±17) 


103 (±0) 


19.2 (±0) 


Bulletrap 


245 


846 (±30) 


143.9 (±3.9) 


17.6 (±0.1) 


Rooifontein 


243 


788 (±9) 


188.9 (±4.7) 


1 7.2 (±0.2) 


Gannabos 


274 


461 (±23) 


204.7 (±8.9) 


18.3 (±0.2) 



Numbers of Aloe dichotoma individuals counted per population and climate characteristics derived for each of the 14 sampled populations in Namibia and South 
Africa. Populations above and below the bold line fall in the summer and winter rainfall zones, respectively. 



conditions at each population from the well established 
and highly resolved (±1 km) interpolated global climate 
dataset, Worldclim (Table 1) [32]. A number of allometric 
measurements were taken for each individual, including 
basal circumference, circumference at first branching 
node, total height and canopy diameter. Trees were then 
subjectively assigned to juvenile, young adult, mature 
adult, senescent or dead age class categories based on allo- 
metric criteria outlined in Table 2. These five categories 
were later reduced to three for ease of analysis and 
interpretability. 

Dead individuals were grouped into windthrown and 
non-windthrown categories based on two main lines of 
evidence. Firstly, windthrown individuals were always 
found lying on their sides, with attached canopies in all 
but the oldest dead skeletons (Figure 1). In addition, no 
alternative explanation could be found for why an indi- 
vidual had fallen over, such as stem or branch damage 
which might suggest an animal had pushed it over, or 
being situated on a very steep and unstable slope where 



toppling over could occur without assistance from the 
wind. Multiple fallen trees orientated in a similar direc- 
tion was not considered an essential criterion for being 
designated as windthrown. This was because an often 
complex topography was likely to cause winds associ- 
ated with discrete storm events to gust in different 
directions, which over time would result in an accu- 
mulation of different fall directions, and windthrown 
trees on even gentle to moderately steep slopes some- 
times rolled downslope, thereby altering their final 
orientation. 

Secondly, windthrown individuals invariably had roots 
attached to the base of the trunk, irrespective of the 
'time since death'. In more recent windthrown examples 
this was often accompanied by soil and rocks still bound 
in the exposed root mass, as well as a depression created 
through the upward leverage of this material. If a tree 
lying on its side had decayed to the extent that the at- 
tached roots were no longer visible, it was excluded 
from the counts. 



Table 2 Age class categories for Aloe dichotoma individuals based on architectural and reproductive characteristics 


Age class 


Original age class 


Plant height (m) 


Canopy 


Dichotomous 


Leaf rosettes 


Reprod. 








diameter (m) 


branching nodes 




mature 


Juvenile 


Juvenile 


Generally <1.5 


ca. 0.60 


None 


1 


N 
















Young adult 


Variable, generally 2-3 


1-1.5 


1-3 


2-8 


Y 


Adult 


Mature adult 


Variable, generally 3-5.5 


2-3.5 


6-10 


>10 


Y 




Senescent 


Generally >4 


2-4 


8-12 


>20 


Y 


Dead 














Dead 


n/a 


n/a 


n/a 


n/a 


n/a 
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Although the possibiUty exists that individuals could 
be blown over after death, and therefore be incorrectly 
classified, this is thought to be a rare occurrence for 
two reasons. Firstly, upon death, individuals soon lose 
water stored in their stems, branches and leaves, mak- 
ing them less top-heavy and therefore prone to top- 
pling over as a result of prior listing or wind gusts. 
In addition, the desiccation of rooting material after 
death seems not to diminish overall rooting strength 
to any great extent. Where individuals have toppled 
over for whatever reason, there is an additional way of 
determining - through inspection of the location of 
the root plate - if the toppling occurred before or after 
death. This is because the root plate is extremely well 
attached to the stem-base in living individuals, but the 
strength of this bond diminishes subsequent to mortal- 
ity. Therefore, a root plate still attached to a long dead 
fallen tree is usually indicative of prior windthrown 
mortality, while an in situ root plate, separated from 
the stem-base, invariably means that toppling occurred 
after death. 

A second reason to question the prevalence of post- 
death windthrow is that leaves and terminal branches are 
lost soon after death, which significantly reduces the drag 
force in the rigid canopy. Dead individuals were placed 
in the non-windthrown category if the dead stems were 



upright and rooted in the ground. The process of non- 
windthrown mortality (i.e. the shift from senescence to 
death) is illustrated in [19]. 

Time since death' 

Estimates of the 'time since death' of both windthrown 
and non-windthrown or standing' dead individuals were 
based on long-term decay rates observed in several repeat 
photographs (M.T. Hoffman, unpublished observations), 
as well as on assessments of decay based on repeated visits 
to individuals known to have died recently. Six categories 
were used to classif)^ 'time since death', using measures of 
structural integrity outlined in Table 3. 

Tree surface area 

Individual tree surface area was calculated using allomet- 
ric measurements and by representing the tree trunk and 
canopy as simplified geometric shapes. The trunk was 
simplified to a trapezium and the canopy to the sector of a 
circle with the triangular base removed (Figure 2). The 
equations used to derive surface areas are presented 
below. Average surface area based on all living adult indi- 
viduals in a population was then graphed to determine 
whether a relationship existed between larger average sur- 
face area in a population and the percentage of wind- 
thrown mortality. 
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Table 3 Criteria used for estimating the number of years since individuals died 



Category 


Years since 
death 


Leaves 


vo DaiK 

remaining 


Branches 


Stem decay 


Growth rings 


vo piani 
remaining 


1 


1-5 


Present 


>95% 


All; intact 


None 


Not visible 


100% 


2 


6-10 


Sometimes; dried out 


50-75% 


All; seldom intact 


None 


Not visible 


90-95% 


3 


11-20 


Absent 


<40% 


All; not intact 


Slight 


Not/barely visible 


80-90% 


4 


21-40 


Absent 


<5% 


>50% lower only 


<25% 


Becoming visible 


-50%, usually trunk 


5 


41-60 


Absent 


Absent 


<10% lowest only 


25-50% or more 


Clearly visible 


-20-40%, trunk only 


6 


61 + 


Absent 


Absent 


Absent 


>80% 


Clearly visible 


<5%, trunk only 



Area of the trunk = (K -h G)T/2 and, 

Area of the canopy = 

{K^-B^) {2HK + S)+ 2B^ {B^ {4H^-K^) + K{4H^K + + 4//^) ArcTan) 



\B-K){B^K) 



2HK + S 



^{b-k)\b^ky 



where, S ^ B^ {K^ -B^ ^ 4H^) , K is the diameter at first 
branching node, G is the basal diameter, T is the height 
to first branching node, B is the canopy diameter, and H 
is canopy height. 




Figure 2 Simplified shape of Aloe dichotoma for determining 
surface area, with dimensions outlined in the text. 



Rainfall, wind-speed and lithology 

A range of climate variables extracted from the interpo- 
lated Worldclim dataset [32] were individually regressed 
against the percentage of windthrown mortality in a popu- 
lation. Variables included mean precipitation of the warm- 
est, coldest, wettest and driest months and quarters, as 
well as precipitation variability. Several measures of wind 
speed were derived from daily wind data for Namibia for 
the longest time period possible, and only included stations 
positioned as close to sampled A. dichotoma populations 
as possible. The most detailed geological and lithological 
maps available were used to describe the rooting substrate 
at each population. 

Modelling procedure 

A logistic regression model with the following form was 
fitted to the available environmental data: 



where p is the probability of death by windthrow, Ci are in- 
dependent climatic variables, dj are independent 'demo- 
graphic' variables (indicators for geology and lithology, 
and basal circumference), and g is an error term. 

To facilitate the interpretation of main and interaction 
effects, all of the continuous independent variables were 
standardized to have zero mean and unit standard devi- 
ation [33]. A stepwise approach was used to find signifi- 
cant main effects, following which all second-order 
interactions between significant main effects were tested 
for significance. Only one of the interactions (between 
mean precipitation of the warmest quarter (PwarmQ) 
and the average of the 5 windiest days per month for the 
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windiest time of day) was significant at a 5% level and 
hence included in the model 

Results 

Percentage of windthrown mortality 

Windthrow was a commonly observed form of mortality 
throughout the latitudinal distribution of A. dichotoma 
(Figure 3), but especially in populations occurring in the 
middle to northern reaches of the summer rainfall zone 
(SRZ) between the Omaruru River and the farm Hau- 
chabfontein (Figure 4). 

Average windthrown mortality across the entire SRZ 
accounted for 31% (n = 11) of total mortality, peaking at 
as much as 72% at Remhoogte. Windthrow incidence 



within the SRZ was much lower at the high elevation 
Brandberg population and south of Hauchabfontein, 
remaining low within the winter rainfall zone (WRZ). 
Here it accounted for only 13% of total mortality, although 
this figure was derived from only three populations. There 
were no windthrown individuals at Rooifontein (Figure 4), 
a population which also registered the lowest number of 
dead trees at thirteen (Figure 3). 

Windthrow and age class 

All age classes had significantly smaller basal circumfer- 
ence measurements compared to the windthrown dead 
category (Figure 5a). However, the mean (M) and 95% 
confidence interval (CI) for the mature adult age class 



22*S- 



23'S- 



24'S 



26'S- 



27-S 



30'S 



31*S- 




Legend 
^ Towns 
Borders 

, Pie charts scaled 
to total no, dead 

I % windthrow 
] % non-windthrow 
Winter/summer 



* rainfall zone 



50 100 



400 Kilometers 



-T — 

12*E 



— 1— 

13*E 



— J— 

14*E 



— 1— 

15*E 



— 1 — 

16'E 



— r — 

17*E 



18'E 19'E 20'E 2rE 
Figure 3 Percentage of windthrown to non-windthrown deaths for each of 14 Aloe dichotoma populations. 
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Figure 4 Percentage windthrown A/oe dichotoma mortality in 
14 populations, arranged from equatorward to poleward range 
extent. Black and white bars fall in the summer and winter rainfall 
zones, respectively. (See Table 1 and Figure 1 for the full names of 
the sampled populations). 



was comparable to that of windthrown dead (M = 197.2, 
CI = [193.4;201.0] and M = 218.2, CI = [210.7;225.7]), 
while it was an order of magnitude smaller for the 
combined juvenile and young adult age class (M = 50.8, 
CI = [47.9;53.6]). In addition, a positive relationship was 
found between the percentage of adult individuals in a 
population and the percentage of windthrown individ- 
uals (Figure 5b), indicating that windthrown mortality is 
more likely to be derived from the larger adult age class. 

Time since death' 

Categorisation of the approximate 'time since death' sug- 
gested that windthrown individuals died slightly more 
recently than non-windthrown or standing' dead, al- 
though both categories indicated that mortality peaked 
some time ago, between 21-40 and 41-60 years before 
present (Figure 6). 
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Figure 5 (a) An illustration of basal circumference differences 
between non-windthrown living and dead individuals and 
windthrown dead and (b) the relationship between the 
percentage of adults in a population and windthrown 
mortality, (a) Bars represent average basal circumferences for age 
classes with standard deviations attached. White bars represent 
non-windthrown living and dead individuals, and black bars 
windthrown dead individuals, (b) Black and white diamonds 
represent summer and winter rainfall zones, respectively. 



70 



Possible drivers of windthrown mortality 

Despite several prior studies linking greater surface 
area to increased windthrow susceptibility, the relation- 
ship between average tree surface area and percentage 
of windthrown individuals in the 14 sampled A. dichot- 
oma populations was poor, with an value of only 
0.07 (Figure 7). 

Wind speeds differed considerably between recording 
stations, especially between coastal and inland stations, 
with Diaz Point (a coastal station) registering the highest 
average and maximum wind speed (Table 4). Conversely, 
the inland station of Mariental had the lowest average 
and maximum wind speed. It was decided to use values 
for the average of the 5 windiest days per month in sub- 
sequent analyses as this was a compromise between an- 
nual average, which would not highlight extreme events, 
and maximum wind speed, which might not always co- 
incide with rainfall events. However, interpretation of 
wind speeds in the context of windthrown death should 




1 2 3 4 5 6 

Index of time since death 

Figure 6 The percentage of windthrown and non-windthrown 
mortality, assigned to an index of 'time since death'. Refer to 
Table 3 for a detailed description of structural criteria used in 
establishing 'time since death' categories. 
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Figure 7 Relationship between average tree surface area of 
mature adult and senescent individuals and percentage of 
windthrown individuals per population. Black and white 
diamonds represent summer and winter rainfall zones, respectively. 



Table 5 Geology (rock type) and llthology (soil type) 
associated with each population 



be made with caution due to short recording periods in 
some cases (e.g. Omaruru), as well as the poor spatial 
network of stations collecting wind speed data and con- 
sequently, the large distances between stations and A. 
dichotoma populations. 

Although there was considerable variability in geology 
and lithology within populations, leptosols predominated at 
most sites (Table 5). Aloe dichotoma populations were also 
strongly associated with rock outcrops in the north and 
with a range of other soil types such as regosols, calcisols, 
arenosols and cambisols, in more poleward populations. 

Measures of wet and dry season rainfall were well corre- 
lated with the percentage of windthrown mortality (posi- 
tively and negatively, respectively), while temperature 
measures were not. Specifically, mean precipitation of the 
warmest quarter (PwarmQ), which is a measure of sum- 
mer precipitation, had a strong positive relationship with 
the percentage of windthrown mortality (R^ = 0.48, p = 
0.01, n = 14; Figure 8a), as did variability in rainfall (R^ = 
0.44, p = 0.02, n = 14; Figure 8b), which is generally greater 
within the SRZ. Since summer rainfall and rainfall vari- 
ability in the WRZ are typically low compared to the SRZ 



Population 

Brandberg 

Omaruru River 

Spitzkoppe 

Tinkas River 

Remhoogte 

Hauchabfontein 

Gorab 

Namtib 

Carolinahof 

Kliphoek 

Grunau 



Rock type 

granite 

complex geology 
granite 

quartzite/granite 

quartzite/limestone 

limestone 

quartzite 

granite 

sandstones/limestones/shales 

sandstone/shale 

dolorite 



Soil type 

leptosols/rock outcrops 
leptosols/rock outcrops 
leptosols/rock outcrops 
leptosols/rock outcrops 
leptosols 
leptosols 

leptosols/regosols 
leptosols/regosols 
calcisols 
leptosols 

leptosols/regosols 



Bulletrap 

Rooifontein 

Gannabos 



gneiss 
gneiss 
shale 



Populations above and below the bold line fall In the 
rainfall zones, respectively. 



leptosols 

arenosols/leptosols 
cambisols 
summer and winter 



(and would therefore influence the slope of the lines and 
strength of the correlations) relationships were assessed 
whilst excluding WRZ populations. This yielded R^ values 
of 0.45 for PwarmQ, and 0.37 for precipitation variability. 
Further exclusion of the anomalous high elevation Brand- 
berg population from the SRZ improved R^ values to 0.81 
and 0.54 for PwarmQ and precipitation variability, 
respectively. 

Logistic regression analysis 

Notwithstanding shortcomings in the data, all available in- 
dependent variables were combined into a logistic regres- 
sion model to assess the relative influence of variables on 
windthrown mortality, as well as their combined predictive 
power (Table 6). Lithological data were chosen over geo- 
logical data due to problems arising from co-correlation. 



Table 4 Namibian wind speed records from stations nearest to Aloe dichotoma populations 



Population 



Nearest station Dist. to nearest No. years Quarter used to Ave. ann. wind Ave. of tiie 5 Max wind 
popn (Icm) on record derive averages speed (m/s) windiest days (m/s) speed (m/s) 



Brandberg/Omaruru River/ Omaruru 
Spitzl<oppe/ 



Tinl<as River/Remlioogte/ 
Haucliabfontein 

Gorab 

Namtib 

Carol inaliof/Kliplioel</ 
Grunau 



Gobabeb 

Mariental 
Diaz point 
Keetmans-lioop 



150/96/ 89 

97/129/ 146 

168 
131 

124/154/151 



3 

23 

12 
26 

26 



NDJ 

JFM 

DJF 
DJF 
DJF 



2.04 

3.37 

1.01 
7.98 
3.91 



5.96 

7.51 

5.18 
18.37 
10.45 



8.06 

9.20 

7.66 
21.02 
12.28 



The quarter from which the average of the 5 windiest days and max wind speed values were derived was based on the months of the warmest quarter (I.e. 
PwarmQ). Values for the average of the 5 windiest days and max wind speed are for the windiest time of day (consistently 20 hOO), with the exception of Diaz 
Point (14 hOO). 
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as well as the assumed greater relevance of the former 
variable. Similarly, precipitation variability was excluded 
because it was highly correlated with PwarmQ, which had 
a stronger tie to the percentage of windthrown mortality. 

Results suggest that our model fitted reasonably, but not 
exceptionally well. A likelihood-ratio test against a null 
model was highly significant = 174.75, p < O.OOl), 



and the pseudo-R^ was 12.0%. Tests for goodness of fit 
and model specification both returned non-significant (i.e. 
favourable) results (Hosmer-Lemeshow goodness-of-fit 
test [34]: xl = 7.25, p = 0.51; link test of squared pre- 
dicted value ;2 =- 0.51, p = 0.41). Using a cut-off of 0.25, 
the model correctly predicted 69% of observed wind- 
thrown dead (model sensitivity) and 67% of observed non- 
windthrown dead (model specificity). The resulting area 
under the ROC curve was 0.73 (compared to 0.5 obtained 
with a null model). 

Odds ratios for main effects can be interpreted as the 
change in the odds of windthrown mortality associated 
with a one standard deviation increase in the independ- 
ent variable. Therefore all variables with an odds ratio 
greater than one (all except regosols) were associated 
with increasing the odds of windthrow mortality by a 
factor of between 1.21 (basal circumference) and 4.86 
(leptosols). Being located on regosols, on the other 
hand, reduced the odds of windthrow mortality by a 
factor of 0.67. 

The significant interaction between PwarmQ and the 
average of the 5 windiest days meant that simultaneous 
increases in these variables made a tree particularly sus- 
ceptible to windthrow. If either PwarmQ or the average 
of the 5 windiest days increased by one standard devi- 
ation while the other was held constant at its mean, the 
odds of windthrown mortality increased 3.90 and 2.59 
times, respectively. However, if both variables were sim- 
ultaneously increased by one standard deviation, the 
odds increased by a factor of 3.90*2.59*3.85 = 38.89. 

Discussion 

Windthrow is a previously unappreciated and important 
source of mortality in A. dichotoma populations. This is 
especially true for populations within the northern half 
of the SRZ, where (with the exclusion of the high eleva- 
tion Brandberg population) windthrow accounted for an 
average of just under half (49%, n=5) of observed mor- 
tality. Even if one accepts a certain degree of error due 
to possible misclassifications of windthrow (as addressed 



Table 6 Logistic model results 



Independent variable 


Odds ratio 


Std. Err. 


z 


P 


95% conf. interval 


Basal circumference 


1.21 


0.08 


2.79 


0.005 


[1.06; 1.39] 


PwarmQ 


3.90 


0.70 


7.63 


<0.001 


[2.75; 5.53] 


Ave. of the 5 windiest days 


2.59 


0.50 


4.92 


<0.001 


[1 .77; 3.79] 


Pwarmq*Ave. of the 5 windiest days 


3.85 


0.90 


5.73 


<0.001 


[2.43; 6.10] 


Leptosols 


4.86 


1.54 


4.99 


<0.001 


[2.61; 9.04] 


Regosols 


0.67 


0.13 


-2.06 


0.039 


[0.46; 0.98] 


Constant 


0.21 


0.05 


-6.39 


<0.001 


[0.13; 0.34] 



The logistic model included continuous and categorical environmental variables, and illustrates the influence of main and interactive effects on 
windthrown mortality. 



Jack et al. BMC Ecology 2014, 14:7 
httpy/www.biomedcentral.com/l 472-6785/1 4/7 



Page 10 of 12 



in the methods section), this still represents a consider- 
able proportion of total mortality. 

With the exception of the south-western corner of 
Namibia, which receives frontal rainfall in winter, the 
SRZ receives its rainfall predominantly through convect- 
ive thunderstorms [20]. These storms originate in the 
north and east as a result of moisture fed into the sys- 
tem by the southerly position of the Intertropical Con- 
vergence Zone during the austral summer season [35]. 
Storm systems typically move southwards and west- 
wards, accompanied by strong winds and heavy localised 
showers [36]. This combination of strong wind gusts 
and associated heavy rain are thought to promote local- 
ised windthrow events through the following mechan- 
ism. After initially strong gusts of wind that herald the 
arrival of the storm, heavy downpours saturate the soil 
and increase its malleability. This is followed by rapid 
absorption and storage of water in the spongy pith of 
the stem, canopy branches and leaves of A. dichotoma 
individuals, making trees (especially large adults) top- 
heavy after heavy precipitation events. The characteris- 
tically rigid canopy structure in A, dichotoma^ coupled 
with shallow rooting depths and/or rooting in unstable 
soil horizons then make individuals prone to wind- 
thrown mortality during strong gusts of wind. Such 
gusts may be associated with the current or - more likely 
- subsequent storm event/s which often arrive before the 
soil has dried sufficiently to lose its malleability. 

Reproductively important mature adults and senescent 
individuals appeared to be particularly susceptible to 
windthrown mortality, judging from the similarity in size 
between windthrown individuals and living mature adult 
or senescent trees, as well as the increased incidence of 
windthrown mortality in populations with proportionally 
more adults. Results from the logistic regression model, 
too, indicated an increase in the odds of windthrown 
mortality as basal circumference increased. 

Notwithstanding the influence of size on windthrow 
susceptibility, there can be little doubt that A, dichotoma's 
more or less rigid architectural structure is also a signifi- 
cant contributor toward the species' predisposition to be- 
ing windthrown. While most trees tend to mitigate against 
the drag forces generated by strong winds through flexing 
in their stems and branches (notwithstanding oscillations 
that may develop in synchrony with wind gusts) and re- 
orientating their leaves to reduce drag [27], this strategy 
for energy dissipation is not possible in A, dichotoma due 
to the relative inflexibility of its stem, branches and leaf 
rosettes. This exposes the species to much greater drag 
forces, especially with respect to larger mature adult and 
senescent individuals with denser canopies, compared to 
smaller unbranched juveniles. 

Compounding this vulnerability is the typically large 
height difference between mature A, dichotoma individuals 



and the surrounding vegetation, as well as the species' 
apparent inability to form closed canopy stands. These 
population level characteristics have been well studied 
in commercial forestry species where tapering vegeta- 
tion heights at the stand-edge and thinning of in-stand 
canopies have been shown to reduce the chances of 
windthrow [1,37]. 

Given the strong indication that windthrown individ- 
uals are derived from the larger mature adult and senes- 
cent age classes [23,24] (Figure 5, this study), the poor 
relationship between the surface area of living mature 
adult individuals and windthrown mortality was surpris- 
ing. This is at least partially attributable to the relatively 
small number of populations sampled, but the greater 
reason probably has to do with micro-topographical and 
-lithological variability within each population, which 
creates pockets sheltered from the wind, and/or areas 
conducive to more secure rooting. It is also conceivable 
that above an as yet undetermined surface area thresh- 
old the likelihood of being windthrown becomes more 
or less equal. In these instances, environmental factors 
(such as those mentioned above) again become the chief 
control in determining whether windthrow occurs, and 
not surface area. 

The vulnerability of adult and senescent A, dichotoma 
individuals within certain populations to windthrown 
mortality is cause for concern, as these age classes, par- 
ticularly the adults, hold a disproportionate amount of 
the reproductive potential in a population. Moreover, 
the characteristically infrequent nature of recruitment, 
estimated at intervals of 10 to 15 years or more in A, 
dichotoma populations growing in more marginal 
habitats [17,19] (M.T. Hoffman & S. Jack, unpublished 
observations), means that these populations are vulner- 
able to windthrow events which might occur on a more 
regular basis. Windthrow can therefore be an important 
influence on population structure, and ultimately, on the 
long term persistence of populations. Smaller, marginal 
populations are at an especially high risk of becoming 
locally extinct as field observations and anecdotal evi- 
dence suggest that it is not unusual for more than a 
dozen individuals to be blown over during a single storm 
event (e.g. pers. comm.: Mr. Robert de Pauw, Kenhardt, 
26/02/2008). 

Of the physical and climatic variables investigated as 
drivers of windthrow, mean precipitation of the warmest 
quarter (PwarmQ) (i.e. the amount of rain which falls in 
summer) and precipitation variability emerged as the two 
single variables best explaining spatial patterns in the per- 
centage of windthrown mortality. In both cases, a partial 
reason for the strong relationship was due to the differ- 
ence in summer rainfall amount and consistency received 
in equatorward summer and poleward winter rainfall 
zones, respectively. But even when WRZ populations were 
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removed from the analysis the relationships were still 
good, and they became excellent when the anomalous 
high elevation Brandberg population was excluded. This 
suggests that the greater percentage of windthrown deaths 
in the Namibian SRZ is attributable, at least in part, to a 
more intense period of summer rainfall, with a likely asso- 
ciated higher frequency and lower return interval of high 
rainfall events. 

The multivariate analysis, which focussed on the SRZ 
due to the higher incidence of windthrown mortality in 
this region, supported this basic conclusion, although 
with some interesting additional findings. When co- 
correlated variables were discarded, the main effects 
contributing to windthrown mortality came chiefly from 
leptosols, PwarmQ and the average of the 5 windiest 
days per month. However, while leptosols had the largest 
odds ratio, this variable had a comparatively poor z- 
score and a high standard error, advising a cautious in- 
terpretation. On the other hand, the higher z-scores and 
lower standard error values associated with PwarmQ 
and the average of the 5 windiest days per month sug- 
gested greater confidence in the odds ratios accompany- 
ing these variables. This indicated that summer rainfall, 
but also accompanying wind speeds, were responsible 
for catalysing windthrow events. Interestingly, despite 
the paucity in wind speed records, when the average of 
the 5 windiest days per month was combined with 
PwarmQ to create a new interaction term, this yielded a 
significant positive effect on windthrown mortality (the 
only significant interaction term). This result, along with 
the positive effect from leptosols and basal circumfer- 
ence, suggested that windthrow events were influenced 
to a large degree by the combined effect of summer 
rainfall intensity and associated high winds, as well as 
local lithology and tree size. While this conclusion 
matches our observations in the field, it is important to 
note that not all of the data are of a good spatial reso- 
lution, and that improvements in this regard would 
greatly aid our confidence in these conclusions. 

This hypothesis for mortality in A. dichotoma is the 
antithesis of previous ideas, which link mortality to 
water balance constraints due to rising temperatures 
[17], or severe and prolonged drought events [19], and 
instead indicates that a high percentage of mortality is 
related to brief periods when there is a surplus of water. 
This study also provides some preliminary evidence that 
windthrown mortality might have a more recent tem- 
poral signature than standing' mortality, but there is lit- 
tle evidence to suggest an increase in more recent 
windthrown mortality associated with increases in ex- 
treme weather events e.g. [21,22]. Instead, the peak in 
windthrown mortality estimated at around 40 years ago 
is likely to be tied to a well documented period of ex- 
tremely wet conditions (and associated storminess) 



during the mid-1970s, which saw widespread flooding 
in southern Africa [20]. Conversely, the peak in non- 
windthrown standing' mortality which occurred as 
much as sixty years ago, may be remnants of particularly 
severe drought conditions which prevailed during the 
early 1930s and 1940s [19,20]. While an interpretation 
of these temporal dynamics is necessarily speculative, an 
increase in the number of sampled populations, more 
temporally consistent sampling, refinements of the 'time 
since death' technique of aging dead skeletons (from, for 
example, additional repeat photographs), and a better 
understanding of the physiological thresholds which re- 
sult in mortality, will help resolve temporal questions 
around mortality, as well as whether windthrown mor- 
tality is increasing or not. 

Conclusions 

The present study has highlighted windthrow as an im- 
portant form of mortality in A, dichotoma, especially 
given that it disproportionately affects reproductively 
mature adult individuals in a species known to recruit 
infrequently. In addition, the study has broadly defined 
the climatic and physical conditions which lead to a 
higher prevalence in windthrown mortality, and illus- 
trated where the above conditions intersect to increase 
the likelihood of windthrow events. However, the avail- 
able data resolution makes it difficult to make spatially 
precise predictions about which populations are likely to 
be most at risk. Moreover, the modelling exercise dem- 
onstrates that threshold conditions required for the oc- 
currence of localised windthrow events necessarily 
involve the combined influence of several variables oper- 
ating in concert, and threshold conditions are likely to 
be dynamic. For example, a decrease in one or more var- 
iables can be compensated for by an increase in others 
with the same net result in terms of probability of wind- 
thrown mortality. Therefore, despite increases in the 
spatial resolution of data, it may be difficult to determine 
where and how threshold conditions will align to cata- 
lyse windthrow. 

Despite these caveats, we can be reasonably sure that 
populations in the equatorward SRZ which receive higher 
amounts and/or frequency of convective rainfall and the 
necessary wind and lithological conditions appear to be at 
a greater risk of windthrow. Populations occurring within 
this latitudinal band may, therefore, benefit from more 
frequent and active monitoring efforts. Smaller, geograph- 
ically isolated and more climatically marginal populations 
dominated by adults should be of the highest monitoring 
priority, as the long term persistence of these populations 
is most seriously threatened. 
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